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Abstract 

We consider a pilot-wave approach for the Dirac theory that was recently proposed 
by Colin and Wiseman. In this approach, the particles perform a zig-zag motion, 
due to stochastic jumps of their velocity. We respectively discuss the one-particle 
O ; theory, the many-particle theory and possible extensions to quantum field theory. 

We also discuss the non-relativistic limit of the one-particle theory. For a single 
particle, the motion is always luminal, a feature that persists in the non-relativistic 
limit. For more than one particle the motion is in general subluminal. 

1 Introduction 

Pilot- wave theories form alternatives to standard quantum theory that are not plagued 
by the conceptual difficulties, such as the measurement problem, that plague the lat- 
ter [1-3]. In the non-relativistic pilot- wave theory of deBroglie and Bohm, systems are 
described by actual point particles which move in physical space under the guidance 
of the system's wave function. This theory has been extended to relativistic quantum 
theory and quantum field theory (for the latter see e.g. [4-6]). In particular, for the 
quantized Dirac field, there exists a pilot-wave approach with positions for particles and 
anti-particles [4] and a Dirac sea approach [7] (see [6] for a comparison). 1 The former is 
stochastic, due to the possibility of particle creation and annihilation, while the latter 
is deterministic. 

Recently, Colin and Wiseman have considered a novel pilot-wave approach for the 
quantized Dirac field [8]. This approach starts from a reformulation of the (massive) 
Dirac theory in terms of a pair of massless Weyl fields of opposite chirality. The Dirac 
mass then yields a coupling of these Weyl fields. As such, they were led to introduce ac- 
tual positions for left- and right-handed Weyl particles, which could be done by adopting 
the particle-anti-particle picture or the Dirac sea picture. In the latter case, the chirality 
of the particles may change stochastically, causing a discontinuity in their velocity. This 
results in a zig-zag motion. In the special case of a single particle, the motion is always 
luminal. In the case of more than one particle, the motion is in general subluminal. 



1 These pilot-wave models introduce positions instead of field configurations. While field configura- 
tions can be straightforwardly introduced for bosons, this seems much harder for fermions [5,6]. 
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Similar zig-zag trajectories have been considered before in different contexts. For 
example, Feynman formulated a path integral approach for the Dirac equation in 2- 
dimensional space-time (called the Feynman checkerboard) , in terms of paths of particles 
that move back and forth at the speed of light [9, pp. 34-36], [10] (see [11] for further 
developments of this model). Recently, also Penrose considered such zig-zag paths as 
Feynman graphs for Dirac fermions [12]. 

Colin and Wiseman were motivated to study these alternative approaches by the 
fact that according to the standard model of particle physics, the fermions are funda- 
mentally massless and acquire an effective mass only through the interaction with the 
Higgs field. Hence, in a fundamental pilot-wave approach, actual positions should be 
introduced for the massless particles and one way to do this is by introducing positions 
for Weyl particles. Not only could this be done by adopting the particle-anti-particle 
picture or the Dirac sea picture, but one can also further choose whether or not to dis- 
tinguish particles of different chiralities. That is, instead of introducing particles that 
are distinguishable by chirality, one may also choose to introduce indistinguishable Weyl 
particles (equivalently one could say that one is introducing positions for massless Dirac 
particles). In the case of the Dirac sea approach, this leads to a deterministic theory 
(which was favored by Colin and Wiseman). 

In the present work, we further explore Colin and Wiseman's zig-zag pilot-wave 
approach to the Dirac theory. In section 2, we start with considering the single particle 
Dirac equation. By decomposing the Dirac spinor into a left- and right-handed chiral 
component, a set of coupled Weyl equations is obtained, for a left- and right-handed 
Weyl spinor. This suggests a pilot-wave model according to which there is an actual 
particle that moves along a continuous trajectory at the speed of light, but with a 
direction of velocity that changes stochastically and either depends on the left or right- 
handed spinor. This model is extended to systems of many particles in section 3. In 
this case, the particles do not necessarily move at the speed of light. In section 4, the 
possible extensions to quantum field theory are discussed. The different possibilities 
arise from the fact that, first of all, as mentioned before, one can adopt the usual 
particle-anti-particle picture or the Dirac sea picture. Secondly, one can regard the 
massive Dirac particles as fundamental or the Weyl particles. Of course, in a more 
fundamental theory, where positions should be introduced for massless particles, one 
should opt for the latter choice. The different possible pilot-wave models in terms 
of Weyl particles are summarized in section 4.4. Then, in section 5 we discuss the 
non-relativistic limit of the single particle model. Finally, in section 6, we consider 
possible pilot-wave models for respectively the Pauli equation and the non-relativistic 
Schrodinger equation that are constructed in a similar way. The latter models should not 
be taken as serious alternatives to the usual approaches, but instead serve to illustrate 
theoretical possibilities. 

2 Single-particle Dirac theory 

Before presenting the zig-zag pilot-wave model, we first consider Bohm's original pilot- 
wave approach [13], which is deterministic. For a single particle, the guiding wave is 
a Dirac spinor ip( x ) = ^( x >*)> which takes values in C 4 , and which satisfies the Dirac 
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equation 2 

ij^D^ip - mip = , (1) 

where = (V, A) is an external electromagnetic potential and = + ieA^ is the 
covariant derivative. In Bohm's approach, the possible world lines x^(s) of the particle, 
parametrized by s, satisfy the guidance equation 

dr^ 

^=fn = iWV - (2) 

where tp = ip^j and jp is the usual Dirac current, which is conserved, i.e. 

d,f D = 0. (3) 

(Note that multiplication of the current by a function a(x) merely changes the parametriza- 
tion of the world-lines and not their image in Minkowski space-time [14] ) . With respect 
to a particular reference frame, the path can be expressed as x(i). The guidance equation 
then reads 

dx j D ip^aip 

^ = VD = S = W (4) 

Since jpjDn > 0, |vd| < 1, so that the velocity never exceeds the speed of light [2]. 
Moreover, for generic wave functions, the probability that the speed of light is reached 
is zero [15]. 

For an ensemble of systems all with the same Dirac spinor ip, Bohm's pilot-wave 
approach reproduces the predictions of standard quantum theory, given that the parti- 
cle distribution is given by ip^ip. This is called the quantum equilibrium distribution. 
The guidance equation implies that if the distribution is given by ip' (x, to)Y>(x, to) at 
a certain time to, it will be given by ip^ (x,t)ip(x,t) at other times t, a property called 
equivariance [3]. This follows from the fact that an arbitrary distribution p(x,t) that is 
transported along the trajectories satisfies the continuity equation 

dtp + V ■ (v D p) = (5) 

and that the density ip'ip satisfies this equation because of (3). 

The zig-zag pilot-wave model forms an alternative approach. It can be obtained as 
follows. 3 Using the chiral projection operators 

Pr = (1 + 75)/2, P l = (1- 75 )/2, (6) 

the Dirac spinor can be decomposed into a right- and left-handed chiral component: 
ipn = PriP, ipL = Pl^P- 4 Using this decomposition, the Dirac equation reduces to 

i^D^ipR = mip L , i^D^ip L = mtp R . (7) 



2 Throughout the paper we use natural units in which h = c = 1. 

3 Colin and Wiseman [8] obtained this model in the context of quantum field theory, adopting the 
Dirac sea picture. We will turn to quantum field theory in section 4. 

4 The decomposition of the spinor corresponds to the fact that the Dirac spinor transforms according 
to a reducible representation of the Lorentz group, given by T>i^' ) © T>( 0, i) [16]. The left- and right- 
handed spinor respectively transform according to the V(h>°) and ©(<>■*) representation. A parity 
transformation transforms one into the other. 
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The Dirac current decomposes as jp = j R + j^, where 

f R = $R-fll>R , 3l = i>Ll^L ■ (8) 

These currents are light-like, i.e. J R jRn = j^JL^ = [8, 16]. They are not conserved, but 
instead satisfy 

d,f R = F, d,f L = -F, (9) 

where 

F = 2mlm (^Vl) • (10) 
The equations (9) can be written as 

OtpR + V • (vrpr) = t LR p L - t RL p R , 

dtPL + V • (v L p L ) = t RL p R - t LR p L , (11) 

where 

PR = J R , PL=J L , v fl = ^ = -*-—, v L = ?r = - T — , (12) 

3r V r Vr 3l W l Wl 

F+ Im (V^Wl) (-F) + Im(^[7Vi?) 

t LR = -o- = 2 ^ \ t ; ' ^ = = 2m ,t ; ' ^ 13 ) 

Jl W l Wl 3r WrVr 

with F+ = max(F,0). 

The equations (11) now suggest a pilot-wave model where the particles move along 
continuous trajectories, with a velocity field that stochastically jumps between or 
vl, with space-time dependent jump rates ti R and t R i to jump respectively from to 
vr and vice versa. 5 Since the currents j R and j£ are light-like, \v R \ = |v^| = 1 and the 
particles always move at the speed of light (a feature that in general does not hold any 
longer in the many-particle case [8]). 

For an ensemble of particles all guided by the same Dirac spinor ip, the velocity 
phase space distribution 

p(x, v, t) = p R (x, t)S(y - v fl (x, t)) + p £ (x, t)5(v - v L (x, t)) (14) 

plays the role of equilibrium distribution. According to this distribution, the probabili- 
ties that a particle is in d 3 x around x, with velocity respectively v R and v £ , are given 
by j R (x,t)d 3 x and j£(x, t)d 3 x. This distribution is equivariant and the dynamics was 
chosen in order to guarantee this (see [4, 17-19] for a general discussion). Namely, for 
an arbitrary distribution 

p(x, v, t) = p R (x, t)6(v - v R (x, t)) + pl(x, i)<5(v - v L (x, t)) , (15) 

the spatial distributions p R and pl satisfy the master equations: 

dtPFL + V • (vrPr) = thRPL ~ tRLPR , 

dtPL + V • (v L p L ) = t RL p R - t LR p L . (16) 



3 Note that Colin and Wiseman also labeled the positions with the chirality. We will not do this here. 
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Because of (11), the distribution given by pr = pR and pi = Pl satisfies this equation. 

The resulting position distribution is pr + pl = ]r + i\ = j% = ip'ip- As such, 
in equilibrium, the predictions agree with those of Bohm's pilot-wave approach for the 
Dirac theory, and hence with the standard quantum predictions. (For non-equilibrium 
distributions, as studied for example in [20], the pilot- wave approaches yield different 
predictions.) 

The theory can be formulated more compactly as follows. By introducing a config- 
uration (x, c) G R 3 x {R, L}, given by a position x and chirality c, the trajectories are 
determined by 

^ =Vc 

and the chirality can jump from c to the opposite chirality, denoted by ire, with jump 
rates i C)7rc . The equilibrium distribution reads p c (x,t). 

A similar approach was proposed before by de Angelis et al. [21] for the case of one 
spatial dimension (who considered it as a generalization of Nelson's stochastic mechan- 
ics [22,23]). In this case, the velocity fields read vr = 1 and vl = —1 (i.e., plus or 
minus the speed of light). The proposed jump rates are actually different from those 
considered here (requiring equivariance of (14) does not uniquely determine the rates). 
The jump rates here are minimal in the sense explained in [4]. 

In the following, it will sometimes be convenient to use an explicit representation for 
the 7-matrices. Here, we choose the Dirac-Pauli representation. In this representation, 
we have, with ip = (jp, x) T , where Ip and % are 2-spinors, that 

1 (f">\. ,i„ = ±( VL 



^=tA4J< (18) 

where ipR = (£> + x)/\/2 and tpi = (<p — X)/V% are Weyl spinors. As such, the wave 
equations (7) become 

la^D^R = mipL , ia^D^ipL = rmp R , (19) 

where = (1, <x) and = (1, —tr). The mass yields a coupling between the left- and 
right-handed spinor. In the case of zero mass, these equations decouple and respectively 
yield the right- and left-handed Weyl equation. 
We further have that 

Or = Pr^PR = (<Pr<PR, <PR<r<PR) > 3l = ¥>i^Vi = (<Pl<Pl, -<p L cr<p L ) , (20) 

F = 2mlm {^ R p L ) , (21) 
so that the velocities and jump rates can respectively be written as 

_, _ 3R _ <Pr°-<pr _, _ h A^PL . . 

Vi? - — - — r , v L - — t , (22) 

3r <p R fR 3 L f L fL 

F+ Im (^Vi) (-F) + Im (^L^i?) 

t L R = — = 2m ^ '— , t RL = = 2m ^ '— . (23) 

3l tp^tpL 3r <Pr<PR 
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3 Many-particle Dirac theory 



The zig-zag pilot-wave model for a single particle can straightforwardly be extended to 
many-particle systems. The wave function ip = ip{x.\, ■ ■ ■ , x n , t) now takes values in the 
n-particle spin space (C 4 )®" and satisfies the n-particle Dirac equation 

n 

id t ip = [~ ia i • ( v i - ieA(xj, t)) + eVfat) + mft] V , (24) 

i=l 

where aj = l<8>---<8>l<8>a:<8>l<8>...l, with a at the ith of the n places in the product, 
and similarly for fa. (For simplicity, we have assumed equal masses and charges for the 
particles.) 

We define the projection operators 

Pc = P C1 ® • • • ® P Cn , (25) 

where c = (ci, . . . , c„), with q = i?, L, which are obtained by taking the product of n 
chiral projections operators. As such, ip = ^2 c i^ c , where ip c = P c ip and the sum is over 
all possible values of the chiralities. Application of the projection operator P c to the 
Dirac equation yields 

n n 

idttpc = ' ( Vi ~ ieA(xi, t)) + eF(xj, t)] ip c + m^2 Mmc , (26) 

i=l i=l 

where 7TjC is obtained from c by changing the ith index q from R to L or vice versa. It 
follows that 

n n 

d t pc + v * • ( v i,cPc) = Yl F w ' ( 27 ) 

i=l i=l 

where 



V'JV'c , v ijC = ^ c ?'^ c , F^ iCjC = 2mlm (^lfaip ni ^\ ■ (28) 

WcWc 



We can now assume tj. particles, with positions Xi, . . . ,x n , whose velocity field is 
given by v c = (vi jC , . . . , v n>c ) for a certain c, and jump rates 

W = ^ (29) 

for the velocity field v c to jump to v WiC . Note that while a jump will only change one of 
the chiralities, it will in general lead to a discontinuity in the velocity of all the particles 
due to entanglement of the wave function (in particular, this will be the case of identical 
particles for which the wave function is assumed to be completely anti-symmetric under 
simultaneous exchange of position and spinor index). The velocity will not exceed the 
speed of light, but unlike the single particle case, it will in general not equal the light 
speed [8]. 

In quantum equilibrium, given by the distribution p c , the position distribution is 
given by YlcPc = V'V) so that this approach yields the same predictions as the usual 
pilot-wave approach to the many-particle Dirac theory (which was presented in [1]). 
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In the Dirac-Pauli representation, we can introduce the many-particle wave function 
ip c which takes values in (C 2 )® n and which is defined as ip C:a = v^^c.a, where a = 
(«i, . . . , a n ), with cti = 1, 2, are the spinor indices. As such the wave equations can be 
written as 

n n 

id t ipc = ^2 i- s ( c i)' lcr i • ( v « ~ ieA(xj, t)) + eV(xi,t)] cpc + m^ Pmc , (30) 
i=l i=l 

where 

s(R) = l, s(L) = -l. (31) 

The equilibrium distribution, the velocity field and transition rates can respectively be 
written as 

t ( , 9 Im i^ C ) ( on, 
Pc = ip' c <Pc , v i)C = s{Ci) — i , t c ^ iC = 2m . (32) 

<fc<Pc Pcfc 

4 Quantum field theory 

As we will see, there are a number of ways to extend this theory to the quantized Dirac 
field. First, one can choose to introduce positions for massive Dirac particles or Weyl 
particles. (Of course, in a more fundamental approach to the standard model where 
particles are considered massless, we should not consider the first option.) Second, one 
can choose to adopt the particle-anti-particle picture or the Dirac sea picture. We will 
briefly explore these different approaches here. We will see that in employing the Dirac 
sea picture there is no difference in the two approaches. On the other hand, there is a 
difference in employing the particle-anti-particle picture. 

4.1 General structure of a pilot-wave model with stochastic jumps 

We start with recalling the general formalism for developing a pilot-wave model for 
quantum field theory which includes stochastic jumps of the configuration [4, 17-19] 
(also called Bell-type quantum field theories). The starting point is the Schrodinger 
equation 

=%(*)) (33) 

for the state vector |^). For simplicity, we will only consider the free Dirac field, for 
which the Hamiltonian is 

H = J d 3 x$(x)(-ia- V + m/3)^(x). (34) 

The Dirac field operators satisfy the anti-commutation relations {Y> a (x), ffi a , (x')} = 
5 aa '5(-K — x') and the other fundamental anti-commutation relations are zero. When 
considering the particle-anti-particle picture, we will always use the normal ordered 
Hamiltonian, denoted by : H :. In this way the vacuum state always has zero energy. 
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The Hamiltonian can also be written as 



H = J d 3 x(ifi[(x)(T ■ V£ L (x) - i^(x)cr • Vfo(x) 

+ m ^(x)^h(x) + ^Jj(x)^l(x) ), (35) 

where the Weyl field operators (p c and (p\, c = R,L, are defined through the operator 
equivalent of (18), i.e. 

(36) 

They satisfy the anti-commutation relations {£> C;a (x), 0\ ,(x')} = <5 cc '<5 aa /(5(x — x'). 

The corresponding pilot-wave model will describe actual point-particles. The number 
of particles may change over time due to particle creation and annihilation. In the case 
there is only one type of particles, the configuration space is given by U n£ N (R 3n x {R, L} n ) . 
An element (x, c) of this configuration space represents a configuration x = (xi, . . . , x n ) 
of positions and a collection of chiralities c = (ci, . . . , c n ). In the case there are differ- 
ent types of particles, the configuration space will be given by the Cartesian product 
of the single particle configuration spaces. In between jumps, the particles will follow 
continuous trajectories determined by one of 2™ velocity fields, labeled by the different 
chiralities. Two types of jumps may occur. First, there may be a jump of just the 
chirality, which leads to a jump of the velocity field of the particles. Second, there may 
be jumps due to particle creation and annihilation. 

The construction of the model requires the choice of a suitable POVM P c {dx). The 
probability that the system with state |vl/(i)) is localized in dx, with chirality c, at time 
t, is given by (^(t)\P c (dx)\^(t)} = p c (x,t)dx. The POVM ^ c P c (dx), where the sum 
is over all possible chiralities, should be an appropriate position POVM, such that the 
usual quantum predictions are obtained in quantum equilibrium. 

The choice of dynamics will now be such that it leaves the distribution p c (x) equiv- 
ariant. It is found by considering the equation 

dt(*{t)\P c {dx)\*(t)) + 2Im{V {t)\HP c (dx)\V(t)) = , (37) 

which follows from the Schrodinger equation. The goal is to bring this equation into the 
form 

d t p c (x,t) + V- [v c (x,t)p c (x,t)] = J dx' [t c t c (x ,x,t)p c ,{x' ,t) - t cc t{x, x' , t)p c (x, t)] , 

d 

(38) 

where the sum is over all chiralities and the integral over all possible position configu- 
rations, and where 

V • [v c (x, t)p c (x, t)} dx = 2Im(*(t)| J ffiP c (da:)|*(t)) , (39) 

Im{y(t)\P c/ (dx')H 2 P c (dx)\y(t))] + 

t cd (x , x , t ) dx' = 2 i j — , (40) 

m)\Pc(dx)Mt)) 
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and H = H\ + H2 corresponds to some natural decomposition of the Hamiltonian. The 
equation (38) then suggest a pilot- wave model with velocity fields v c (x, t) and jump rates 
t cc '(x, x', t)dx' for a configuration to jump from x to an element dx' around x' and from 
c to d . 

We may obtain a pilot-wave model with no distinction between different chiralities, 
by starting from the position POVM ^2 c P c (dx), where the sum is over all possible 
chiralities. (See [4, 19, 24] for a general discussion on how to construct a pilot-wave 
approach in terms of identical particles.) 



4.2 Massive Dirac particles 
4.2.1 Particles and anti-particles 

Consider the usual plane waves 

n s (p)e i P- x , v s (p)e~^, (41) 

which are eigenstates of the Dirac Hamiltonian — ict-V+m/3 with eigenvalues respectively 
E p and —E p , where E p = \/|p| 2 + m 2 , and with momentum respectively p and — p (see 
e.g. [25]). The label s = 1,2 denotes the helicity. The Dirac field operator can be 
expanded in terms of these: 



d 3 p 



m 



6 i (p)« a (p)e i P- x + ^(pK(p)e- i P- : 



(42) 



The operators 6 s (p) and d s (p) respectively annihilate a particle and an anti-particle 
with momentum p and helicity s. 

Writing V( x ) = ^( x ) + <^( x )) where b and d) are respectively the particle and the 
anti-particle part of the field operator ip, the Hamiltonian reduces to: 

H = ^2 [ d 3 x [ftt ( x ) + d Q (x)] (-ia • V + m/3) a/3 [^(x) + dj(x 

= ^2 d 3 x %(x)(-ia- V + m/3) a ^(x) + cT Q (x)(-ia- V + m/3) aj8 dj(x) 

a,/9 J 



(43) 



where 



a,/3 



Hd,i = d 3 x fa (x)(-ia a/3 . V)6 /3 (x) + cT a (x)(-k* a/3 • V)dJ(x) 



-H"D,2 = / ^ X ^U x )/WfyK X ) + ^(x^^C^x) 



(44) 



In the second line of (43), we have used the fact that the cross-terms containing par- 
ticle and anti-particle operators vanish (which essentially follows from the property 
ut(p)v s >(— p) = [25]). This leads to a decoupling of the particle and anti-particle part 
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in the Hamiltonian. In the Schrodinger equation, we use the normal ordered Hamiltonian 
: H :. 

Consider now the states 
\x,a;n,n) D = -^=^1^) . . . \ n (pc n )d\ n+1 (x n+ i) . . . d\ (x n+n )\0) D , (45) 

where n, n G N, x = (xi, . . . , x n+s ), a = (a\, . . . , a n+s ), with Qj = 1, . . . , 4, are the 
spinor indices, and \Q)d is the state which contains no particles or anti-particles, i.e., 
6 Q (x)|0)z) = <i Q (x)|0)£> = for all x and spinor indices a. These states correspond to n 
particles and n anti-particles, and span the Hilbert space. They are anti-symmetric un- 
der simultaneous exchange of positions and spin indices corresponding either to particles 
or anti-particles. In [4], a pilot- wave model was constructed starting from the position 
POVM Pp' n \dx) = J2 a \x,a;n,n)DD(x,a;n,n\dx. In this model, there are positions 
for particles and anti-particles. In the absence of interactions, there are no jumps. In 
the case of for example electromagnetic interaction, there are jumps corresponding to 
particle creation and annihilation. 

Here we consider a different POVM. Let us therefore first introduce the states 

\x, c, a; n, n) D = ^ (Pc) a ,p K /3; n, ™)d , (46) 

P 

where c = (c±, . . . , c n+ n) and P c is the projection operator defined in (25). For a state 
\^(t)), the Schrodinger equation (33) (with a normal ordered Hamiltonian) implies that 
the expansion coefficients ^c^\x,t) = £i(x,c,a;n,n\^(t)} satisfy the many-particle 
wave equation (26) for vanishing electromagnetic potentials. 

We can now define the POVM P^f \dx) = £ Q | x, c, a; n, u)dd{x, c, a; n, n\dx. In 
the corresponding pilot-wave model, the velocity field corresponding to : Hp i : is the 
one given in (28), evaluated for the wave function ipi^^fa, t), and the jump rates, which 
are derived from : Hb,2 are 

(nn) \imU^\x,m4 n ' n \x,t))} + 

t ( ":r(x,x',t) = 5^. cd 5(x - x')2m± v , _., '-L- . (47) 

Because of the delta-function S(x — x'), there are no jumps of the position configuration. 
In particular, there is no particle creation or annihilation (even though the wave function 
need not be a particle number eigenstate). There are only jumps of the velocity with 
rates that are the same as those introduced for the many-particle Dirac theory. As 
such, this pilot-wave model reduces to the one for the many-particle Dirac theory. This 
is of course an artifact of the free theory. In the case of for example electromagnetic 
interaction, particle-anti-particle pair creation will be possible, which will be reflected 
in the jump rates. 

In this pilot-wave approach, particles and anti-particles were introduced as different 
particle species. One could also construct a model in which they are indistinguish- 
able. Such a model could be obtained by starting from the alternative POVM which is 
obtained from Pc U,n \dx) by summing over n and n, keeping n + n fixed [4, 19, 24]. 
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Note that in order to have a pilot-wave model which does not distinguish particles 
of different chirality, one could start from the alternative POVM obtained by sum- 
ming Pp'™\dx) over all the possible chiralities c. The resulting position POVM is just 

Pp' n \dx), which was considered in [4]. 



4.2.2 Dirac sea 

Instead of using the particle-anti-particle picture, one can also adopt the Dirac sea 
picture. The starting point is the field expansion 



y>(x) 




b s (p)u s ( P )jv* + b s ( P )v s (- P yv- 



(48) 



where b s (p) = c/J(— p) is the annihilation operator for a particle of negative energy E p 
and momentum p. 
The states 

\x, a) D = -LA (xi) ■ ■ ■ $ an (x n )|0) D (49) 



can be introduced, where the state \Q)d is now the state that does not contain particles 
of positive or negative energy, i.e., V'aMlCOD = for all x and spinor indices a. They 
are anti-symmetric under simultaneous exchange of positions and spin indices and span 
the Hilbert space. The state \0)d can be obtained by applying all the creation operators 
of negative energy particles to \0)d- 

In [7], the position POVM Po(dx) = ^ a \x, cx)dd{x, a\dx was considered. In the 
corresponding pilot-wave model, there is a fixed number of particles which move de- 
terministically, even in the case of electromagnetic or other fermion number preserving 
types of interaction. 

Consider now the states 

\x,c,a) D = (P c ) Q)/9 \x,/3) D . (50) 

The expansion coefficients ip c ,a(x, i) = d{x, c, a\^(t)) of a state !*(£)) satisfy the many- 
particle wave equation (26) as a consequence of the Schrodinger equation (33). Starting 
from the position POVM P£) tC (dx) = ^ Q \x,c,a)oD{x,c,a\dx, we obtain the following 
pilot-wave model. The velocity field is determined by Hd,i and corresponds to the one 
given in (28) for the wave function ip c ,a(x,t) and the jump rates derived from Hd,2 are 
of the same form as in (47). So again, the pilot- wave model is formally similar to the 
one for the many-particle theory. This is still the case when electromagnetic interactions 
are considered, unlike the zig-zag model obtained in the particle-anti-particle picture. 
Another difference with the latter model is the number of particles. For example, the 
case where there are no particles or anti-particles corresponds to a filled Dirac sea and 
hence to an infinite number of particles in the present model (employing regulators the 
actual number of particles could be made finite [7]). 



11 



4.3 Weyl particles 



4.3.1 Particles and anti-particles 

Instead of expanding the field operators in terms of eigenstates of the Dirac Hamiltonian, 
they can also be expanded in terms of eigenstates of the right- and left-chiral Weyl 
Hamiltonians =Fier • V [8]. By introducing the 2-component spinors iOj(p), i = 1,2, with 



<7-p 



w i (p) = (-l) 1 + i w i (p), 



we have that the plane waves 



wi(p)e ipX , ^ 2 (- P )e- ip - > 



(51) 



(52) 



are eigenstates of the right-chiral Weyl Hamiltonian with eigenvalues respectively |p| and 
— |p|. Since they have momentum respectively p and — p, the relations (51) imply that 
they have respectively right- and left-handed helicity, i.e. their momentum direction 
is respectively aligned and anti-aligned with their spin. (More generally, for positive 
energy plane waves, the chirality equals the helicity. For negative energy plane waves, 
they are opposite.) Similarly, we have eigenstates of the left-chiral Weyl Hamiltonian. 
In terms of the plane wave expansion, the field operators read 



d 3 p 



(x) = / 



b R (p) Wl (p)e ipx + d t L (p) W2 (-p)e- ip - x 
6 L (p)u; 2 (p)e ip - x + 4(p) Wl (-p)e- ipx 



(53) 



where 6 c (p) and d c (p) respectively annihilate Weyl particles and anti-particles with 
momentum p and chirality c. They satisfy the usual anti-commutation relations. 
We can also write 



<Pr(x) = b R (x) + d\(x) , ^l(x) = b L (x) + c^(x) . 



51 



(54) 



Using these expressions, the Hamiltonian H can be written as Hw,i + Hw,2, where H\y,2 
corresponds to the term in H which depends explicitly on the mass. We have that 



H 



WA 



b L,J X )( [(T aP ■ V ) fo L,/3( x ) + d R,a (x) (i<X Q(3 • V)^(x) 
-^,a( X )( icr a/3 • V)6 R , /3 (x) - d L , a (x)(iff a p ■ V)d£ ij9 (x)J , 

H w ,2 = mJ2 / d z x (&I !Q (x)^ )a (x) +V Ra d f Ra (x) + d LiQ (x)fe LiQ (x) + dk,a(x)6fl, a (x)) , 

(55) 

where we have used certain properties of the spinors iOi(p) given in [8]. Note that 
H\y,i / Ho,i {Hd,2 is not just the term in H which explicitly depends on the mass). We 
further apply a normal ordering to the Hamiltonian. Note that the Hamiltonian : H\y,2 ■ 
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does not commute with the total particle number. Pair creation and annihilation will 
be possible of a particle and an anti-particle with the same chirality. 
Let us now introduce the states 



x, c, a; n, n)w = 



Vn\n\ 



F^ 1 ,a 1 (xi)...6t n)Qn (x n )4 



'c n +i,a n +i 



(x n+ i)...4. 



(yLn + n)\0} W , (56) 



where a = (cci, . . . , a n +n), with on = 1,2, and \0)w is the state that contains no left 
or right-handed chiral particles or anti-particles. Note that \0)w 7^ |0)_d; \0)d is an 
eigenstate of the Hamiltonian : H : (with eigenvalue zero), whereas \Q)w is not [8]. The 
overlaps w{x, c, a; n, re|\I/(i)) = (^"^(x, t) do not satisfy the many-particle equation 
(30), since the Hamiltonian : H : does not commute with the total particle number 
operator. 

A possible choice of POVM is P$f\dx) = £ Q | x, c, a; n, n)ww{x, c, a; n, n\dx. The 

velocity field corresponding to : Hw,i '■ is the one given in (32). The rates to jump 
from (x,c), with n particles and n anti-particles, to (x',c'), with n' particles and n' 
anti-particles, are given by 



We will not present the explicit expression for these jump rates. There will not be 
jumps solely of chirality. But there may be jumps corresponding to pair creation or 
annihilation of a particle and an anti-particle of the same chirality. 

4.3.2 Dirac sea 



The operators 6r(p) and &l(p) annihilate a particle with momentum p, with negative 
energy — |p| and with respectively right- and left-handed chirality (and left- and right- 
handed helicity). 

The Hilbert space is spanned by the states 




(57) 




x, c, a)w = 



1 



£Ji,ai( x l) ■ ■ ■ <Pl„,a n ( x n)\0)w , 



(59) 



\fn\ 
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where \0)w is the state that contains no Weyl particles of positive or negative energy, 
i.e., (p c ,a(x)\0)w = 0, for all x, c = R, L and a« = 1,2. Due to (36), \0)w = |0)r> and 
\x, c, a)\y = y/2™\x, c, a)o for a = 1,2. The overlaps tp c ^ a {x,t) = w(x,c,a\ i S(t)} satisfy 
the many-particle wave equation (30) . 

Choosing the POVM Pw c (dx) = ^2 a \x, c, a)ww( x i c > ct\dx = Pd c (cIx), the pilot- 
wave model is the same as the one obtained in section 4.2.2 starting from the Dirac sea 
picture for massive Dirac particles. This equivalence still holds in the case of interactions. 

Starting from the position POVM ~51 c Pw,c{dx) = Po(dx) the deterministic pilot- 
wave model of [7] is obtained. So in this approach, particles are not distinguished by 
chirality. One could also say that in this model positions are introduced for massless 
Dirac particles [8] (because the operators 6 c (p) can also be interpreted as the annihilation 
operator of a positive energy massless Dirac particle with helicity c and similarly for the 

operators 6 c (p)). 
4.4 Summary 

We have discussed a number of possible pilot-wave approaches for the quantized Dirac 
field. In the context of the standard model, where the particles are fundamentally 
massless, one should of course introduce positions for the massless particles. We have 
discussed four ways of doing this. One can choose to adopt the particle-anti-particle 
picture or the Dirac sea picture. One can further choose whether or not to distinguish 
particles of different chirality. In the particle-anti-particle picture, there is the possi- 
bility of pair creation, even in the free case, whether or not one chooses to distinguish 
particles of different chirality. In the Dirac sea picture, there are no jumps of position. 
If one introduces particles with different chirality then there might be jumps of chirality, 
otherwise the model is deterministic (even in the case of fermion number preserving 
interactions) . 

5 Non-relativistic limit 

In this section, we derive the non-relativistic limit of the zig-zag pilot-wave model for 
the case of a single particle. 

5.1 Non-relativistic limit of the Dirac equation 

We first consider the non-relativistic limit of the Dirac equation, following the usual 
account as found in e.g. [26,27]. Starting from the Dirac-Pauli representation and the 
ansatz 

*= (!) 

and assuming \idtx ~ e ^x\ ^ m \x\, X can De expanded as 

X = Xi + X3 + • • • , (61) 

where 

^ = - ^ ■ D ^< * 3 = - • D)3 ^ - £^ • ^ 
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with D = V — ieA. E = — <9fA — VV and B = V x A are respectively the electric 
and magnetic field. As such, \ is the small component relative to p, with \xi\ and |%3 1 
respectively of the order (p/m)\ip\ and (p 3 /m 3 )\p\. The dots in (61) denote the higher 
order terms. 

To lowest order, elimination of x m favor of p yields the Pauli equation: 



1 e 
\d t p = -—D 2 p - - — B • ap + eVip . 
2m 2m 



(63) 



The next order term in the Hamiltonian is of the order p 4 /m 3 and only arises when \3 
is taken into account. Those higher order terms will not be considered. 
The expansion (61) can be used to expand the current 



Ail = AV 



(64) 



where c = R,L and the terms are of the order p % jm % . This yields 

•0 1 t 

Jc,o = 2 < PV > > 



7° 



3c,2 



8m 2 



{a ■ T)p) ] (a ■ Dp) 



1 



8m 2 L 



(Dp) j • (Dp - ier x Dp) 



(65) 



and 



Jc,0 = s(c)-^<Tip, 



1 



jc,2 = (tr ■ Dp)' a (a ■ Dp) 



s(c) 



8m 2 
1 

8m 2 



(cr • Bp) ] Dp + (D^) f cr Tip - (Dp)'' • Der^ - i (Dip)* x D ip] , (66) 



with s(c) as defined in (31). Similarly, F, which is given in (21), can be expanded into 
terms Fj which are of the order p l jm l ~ x . The lowest order terms read 



^0 
Fi 
F 2 
F3 



0, 

0, 



Am 2 



Re 



p^D • (cr • D)D(/? 



+ -BV 



+ — E- ^ . (67) 
2m 



Assuming that p satisfies the Pauli equation (63), we have the following identities: 

V-j C)0 = s(c)F 1 , (68) 
%' c °o + V-j c ,i = 0, (69) 
% c a + V-j C)2 = s(c)F 3 . (70) 
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This can be verified by direct calculation. Alternatively, it can be seen by considering 
the expansion of the identity d^jc = s(c)F using (64)-(67). The equations (68)-(70) 
correspond to the terms in this expansion of the order p l+l /m l , with i = 0, 1, 2. Namely, 
the terms V ■ j Cj j are of order p t+l jm % and, assuming the Pauli equation to evaluate the 
time derivatives, the terms dtj^ are of order p l+2 /m t+l . Higher order corrections to 
the Pauli equation, which are of the order p 4 /m 3 , only show up in higher order terms 
in dujc = s(c)F, which are not considered here. 
We also have the identities: 

0c,o) 2 = Lkol 2 , Jcfijc,! = jc,0 • Jc,l , {jc,l) 2 + 2j° j° 2 = |j c ,l| 2 + 2j Ci0 ■ j c ,2 • (71) 

As before, they can be found by considering the expansion of the identities jcjc^ = 
using (64)-(67). 



5.2 Non-relativistic limit of Bohm's approach 

Let us now consider the non-relativistic limit of Bohm's approach to the Dirac theory 
(which was derived before in [1]). In this approach, the velocity field of the particle is 
given by vp, = in/ 3%-, with the Dirac current defined in (2). Since jp = jfi + j^, we 
have from (65) and (66) that 

Ao = ¥>V, Ai = ° ( 72 ) 

and 

Jd,o = 0, j A i = — Im (</>W) + ^ V x (<pW) , JD,2 = 0. (73) 

Hence, up to order p 2 /m 2 , the Dirac current }d is given by This forms a current 
for the Pauli theory and we will denote it by jp. Namely, assuming the Pauli equation, 
it satisfies the conservation equation 

$(¥>V) + V-j P = (74) 

(which is just a rewriting of (69)). 

Up to order p 2 /m 2 , the Dirac velocity field reduces to the Pauli velocity field 

vp = -^f - = — ^Im U-Dip) + — ^ V x Utrtp) . (75) 

The Pauli velocity field can be used to yield a pilot- wave approach to the Pauli equation. 
Namely, in virtue of the continuity equation (74), the corresponding particle dynamics 
guarantees the equivariance of the Pauli distribution 93^. 



5.3 Zig-zag pilot- wave approaches for the Pauli theory 

Before passing to the non-relativistic limit of the zig-zag pilot-wave approach for the 
Dirac theory, it is worth considering the following two pilot-wave approaches for the 
Pauli theory. They are suggested by considering the equations (9) respectively up to 
order p 3 /m 2 and p 2 /m (which correspond to combinations of the identities (68)-(70)). 
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In the next section, we will compare these approaches to the non-relativistic limit of the 
zig-zag pilot-wave approach for the Dirac theory. 
In the first combination of (68)-(70) yields 



dt (j c ° + i c °i) + V • (j c , + j c ,i + j Cl2 ) = s(c)(Fi + F 3 ) 



(76) 



where c = R, L and where the Pauli equation is assumed. This suggests a pilot-wave 
approach with velocity fields 

~ _ jc,0 + jc,l + jc,2 



and jump rates 



tc,TVC 



7° 4- 7° 

Jcfl + Jc,l 



-s(c)(Fi + F 3 )]- 

Jcfl ^ Jc,l 



(77) 



(78) 



for the velocity field to respectively jump from v c to v nc . The equivariant equilib- 
rium distribution is given by + (which is positive under the assumption that 
ico ^ | Jc 1 1 ) ■ The resulting position density is j^, + j° L = <^Vj so that this approach 
reproduces the usual predictions of the Pauli theory. 

In the second combination of (68) and (69) yields 



%c°o + V- {ic,o+}c,i) = s(c)F 1 
This suggests a pilot- wave approach with velocity fields 

- _ jcfl + jc,l ./.a^' "^ , JP 



7° 



*(c) 



+ 



s(c)s + Vp 



and jump rates 



[-s(c)^i] + _ [- S (c)V-(^W)]' 



7° 
Jc,0 



(79) 



(80) 



(81) 



where the vector s = tp^cnp/ip^ip is the normalized spin polarization vector and vp 
is the usual Pauli velocity field. The equivariant equilibrium distribution is given by 
Jcfl = *P and again yields ip'ip as position density. So also this approach reproduces 
the usual predictions of the Pauli theory. 

5.4 Non-relativistic limit of the zig-zag approach for the Dirac theory 

Turning to the zig-zag approach for the Dirac theory, we have that, up to terms of the 
order p 2 /m 2 , the velocity fields read 



1 / ic l 

v c = — jcfl + jc,l + jc,2 rff" (jc,0 + jc,l) + 

Jcfl \ Jcfl 



and up to terms of the order p 3 /m 2 , the jump rates read 



4 1 

. J c,0 



7° 
7° 

Jc,0 



Jc,0 +0 



(82) 



* c ,™= — ( [s(c)(F 1 + F 3 )} + + 

Jcfl 




Jc,l J_ Jc,2 
7° 



(- S ( C )iq)+ +0 ^ 



(83) 
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Let us now compare this to the first model for the Pauli theory considered in the 
previous section. When expanding v c to the same order, we find the same expression, up 
to the term containing j° 2 (which is of the order p 2 /m 2 ). Similarly, when expanding i C 7rc 
to the order p 3 /m 2 , it is equal to the above expression, up to the the term containing j° 2 . 
So by considering the velocity fields and the jump rates up to the present order, we do 
not obtain exactly a pilot-wave approach to the Pauli equation. Namely, the dynamics 
does not ensure the equivariance of the distributions j® + j®^ exactly, but only up 
to terms of order p 2 /m. As such, the position distribution ip^ip is only approximately 
preserved. 

Note that the extra term in the expansion of v c compared to v c follows from the fact 
that |v c | = 1. The expression in the right hand side of (82) has norm one, up to order 
p 2 /m 2 (which can be verified explicitly using the identities (71)), whereas v c does not. 

A similar conclusion holds if we consider the non-relativistic limit up to a lower 
order. Namely, up to terms of the order p/m, we have 

1 /. x . Jc°i. \ Mn fp 2 

v c - t„- Jc,o + Jc,i - ^J c ,o + U I —. - 2 

Jcfl \ Jcfl / ^ 



v c + v A + O 



p 2 
m 2 



•- •• s • v.. v : c > ( ) . (81.) 



where (80) was used and 



Jc,iJ c ,o 1 Im (V? f c- • D99) rfaip 1 Im (y?t<r • D99) 

v A = -— — = 1 r~ = 1 s • ( 85 ) 



The jump rates reduce to 



_ (-s(c)F 1 y ( f c A (jfl 



3c,0 \ Jc,0 



^fi-»(c)l Im(yt r Dy) 




(86) 



The additional term va compared to v c is of order p/m and ensures that s(c)s+vp + 
va has norm one, up to order p/m (which can be verified explicitly using (71)). Note 
that we can also write = — (vp • s) s (using (71)). This implies that the additional 
terms in the velocity fields and the jump rates only vanish when the usual Pauli velocity 
field is orthogonal to the spin vector. 

Note that the dominating term in the velocity field v c is j c ,o/ic = s ( c ) s > which has 
norm one. So up to lowest order of approximation, particles move back and forth at the 
speed of light, along the direction of the spin vector. 

Again we do not obtain exactly a pilot-wave approach to the Pauli equation. Namely 
the dynamics does not exactly preserve the position distribution ip^<p. 
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5.5 Spin eigenstate 

For the special case of vanishing electromagnetic potentials and a Pauli spinor which is 
a spin eigenstate, i.e., v?(x, t) = ^/>(x, where ip is a scalar function and = 1, the 
Pauli equation reduces to the non-relativistic Schrodinger equation for ip. 

Writing ip = \^\e lS and using s = £t<r£, the lowest order terms in the expansion of 
the current and the quantity F obtain the following form: 

ic o = ^H 2 , (87) 

& = s(c)^M 2 8-V5, (88) 

Jc°2 = ^ [IW| 2 + is • W x VV)] , (89) 

Jco = S (c)|s|^| 2 , (90) 

jc ' 1 = i^ |2V5+ 4^ Vx(s| ^ ' (91) 
j c2 = s(c)^ t s • W*VV> + V^*s • W - s|VV>| 2 - iVV>* x W] , (92) 

F = , Fi = • V|V| 2 , F 2 = , F 3 = ^ Re (V>*s • VV 2 V>) . (93) 

These expressions can be used to obtain the form of the zig-zag pilot-wave models for 
the Pauli theory considered in section 5.3 and the non-relativistic limit of the model for 
the Dirac theory. 

Let us first consider the zig-zag model for the Pauli theory defined by (80) and (81). 
The velocity fields reduce to 

v c = s(c)s + — VS+-^-V x (sln|V| 2 ) , (94) 

and the jump rates to 

t c ,*c= [s(c)s- Vln|V| 2 ] + . (95) 

So the velocity field is given by the usual velocity for the non-relativistic Schrodinger 
theory, namely VS/m, plus a curl term (which needs to be added to the usual velocity 
field when one starts from the usual pilot-wave approach to the Dirac theory [28]), and 
plus or minus the spin vector. The jump rates are smaller when the spatial variation of 
density |^| 2 is slower. 

In the case of the zig-zag pilot-wave model for the Dirac theory, the velocity fields 
reduce to 



and 



v c = v c --s(s.VS) + o(4) (90) 



and the jump rates to 



{] ~ s{c) i s ' vs ) +0 d ] (!)7) 
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The term in the velocity field v c that is additional to v c corresponds to minus the 
component of the usual de Broglie-Bohm velocity field V5/m for the non-relativistic 
Schrodinger theory along the direction of the spin vector. 

Similarly as in the case of a general Pauli spinor, the expressions (94) and (95) 
determine a pilot-wave model for the non-relativistic Schrodinger theory because they 
preserve the density |V>| 2 , while the expressions (96) and (97) do not. 

6 Variations on a theme 

We have explored the zig-zag pilot-wave model for the Dirac theory. This model forms 
a stochastic alternative to Bohm's original deterministic approach. For a single particle, 
the stochasticity arises from the fact that the velocity field might jump between two 
alternatives, respectively determined by the left- and right-handed chiral component of 
the Dirac spinor. It is clear that still other pilot-wave models could be obtained by 
starting from a different decomposition of the Dirac spinor. For example, one could 
decompose it into two Majorana spinors, or one could decompose it into the two spin 
components along a certain direction. Since such approaches seem rather unnatural, we 
will not pursue them further. We will just illustrate certain aspects of such alternative 
models in the context of the Pauli theory and the non-relativistic Schrodinger theory. 
Again, there seems to be no reason to consider these models as serious alternatives to 
the usual models. 

6.1 An alternative pilot-wave model for the Pauli theory 

Reconsidering the Pauli equation (63), a different pilot-wave model than the usual one 
can be obtained as follows. Writing 

with (pi and ip2 respectively the spin-up and spin-down component along the z-axis, the 
Pauli equation implies the following identities: 

dt\<Pi\ 2 + Vji = /, 9^ 2 | 2 + V-j 2 = -I, (99) 

where 

j, = -InuV*Dc^) , * = 1,2, (100) 
m 

and 

/ = — Im ((Si + i5 2 )v?2>i) • (101) 
m 

This suggests a pilot-wave model, in which a point-particle moves along continuous 
trajectories, with a velocity field that jumps between vi = ji/|^i| 2 and v 2 = j 2 /|<£> 2 | 2 , 
with jump rates i 2 i = I + /l^] 2 and fi 2 = (— I) + /\ifi\ 2 to jump from v 2 to Vi and vice 
versa. As such, the particle is either guided by the spin-up or spin-down component. As 
a distribution on velocity phase space, the equivariant equilibrium distribution is given 
by 

p(x,v,t) = |( / 9 1 (x,t)| 2 <5(v-v 1 (x,t)) + | V 9 2 (x,t)| 2 5(v-v 2 (x,t)). (102) 
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The corresponding position distribution is tp^ip = \(fi\ 2 + |^2| 2 , so that the model is 
empirically equivalent to the usual model for the Pauli theory. 

Note that in the special case of vanishing electromagnetic potentials and a spin 
eigenstate <p(x, t) = ip(x, this model reduces to the usual non-relativistic pilot-wave 
theory of deBroglie and Bohm for spinless particles. Namely, in that case, the jump 
rates are zero and the velocity field is given by VS/m, where ip = l^le 15 . 

Since we started from a decomposition of ip into its spin components along the z- 
axis, the obtained pilot-wave model is not rotationally invariant. In [29] a similar theory 
was proposed in the context of Nelson's stochastic theory. There it was noted that, in 
the case of a homogeneous magnetic field, rotational invariance could be achieved by 
decomposing the spinor into its spin components along the direction of the magnetic 
field. For a non-homogeneous magnetic field one could consider a decomposition of 
tp> that varies spatially as well as temporally. However, the magnetic field could be 
zero in certain regions of space, leaving such a decomposition ill-defined. A better 
possibility would be to consider a decomposition along the direction of the spin vector 
s = ip^ (Tip I ip^ ip which is non-zero when ip is non-zero. However, since such a construction 
appears rather unnatural and complicated, we will not pursue it further. 

6.2 An alternative pilot-wave model for the non-relativistic Schrodinger 
theory 

Consider now the non-relativistic Schrodinger equation 

id t ^ = — — D 2 Y> + Vij> . (103) 
2m 

We can write ip = tp± + iip2, where ip\ and ip2 are real. This leads to 

^ 2 + V-j 1 = /, d^l + V -h = ~I, (104) 

where , 

ji = —tpi (Vip 2 - eAtpi) , j 2 = tp2 (VV>i + eAtp 2 ) , (105) 

m m 

and 

1 e 
I = —Vtp! ■ VV> 2 + — |A|ViV> 2 + 2FViV>2 • (106) 
m m 

This suggests a pilot-wave model where the velocity jumps between vi = ji/V^i an d 
v 2 = J2/V'2' w ith jump rates t 2 i = I + anc l = ( — -0 + /^i- Note that in this case 
the velocity fields depend on both components of the wave function. The equilibrium 
distribution is given by 

p(x, v, t) = Vi (x, t) 2 5(v - Vl (x, t)) + V 2 (x, t) 2 5(v - v 2 (x, t)) , (107) 

with corresponding position distribution \tp\ 2 . 

The obtained model is not gauge invariant. Namely, a gauge transformation ip — > 
e ip, A —7- A — V# will yield different trajectories. Nevertheless, in equilibrium, the 
violation of gauge invariance can not be detected. In the usual pilot-wave approach of 
deBroglie and Bohm, particles move with a velocity Im(Vip/t())/m = (ji +h)/\' l P\ 2 i 
which is invariant under gauge transformations. 
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The decomposition of the Dirac spinor into two Majorana spinors (which corresponds 
to a splitting into real and imaginary part in the Majorana representation of the Dirac 
matrices) is similarly not invariant under £7(1) gauge transformations. Hence, a pilot- 
wave model in the spirit of the zig-zag model, starting from such a decomposition, would 
also not be gauge invariant. 
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